Study of the ground state properties of LiHo x Yi_ x F 4 using /iSR 
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LiHo x Yi„ x F4 is an insulating system where the magnetic Ho 3+ ions have an Ising character, 
and interact mainly through magnetic dipolar fields. We used the muon spin relaxation technique 
to study the nature of the ground state for samples with x=0.25, 0.f2, 0.08, 0.045 and 0.0f8. 
In contrast with some previous works, we have not found any signature of canonical spin glass 
behavior down to ~I5mK. Instead, below ~300mK we observed dynamic magnetism characterized 
by a single correlation time with a temperature independent fluctuation rate. We observed that this 
low temperature fluctuation rate increases with x up to 0.08, above which it levels off. The 300mK 
energy scale corresponds to the Ho 3+ hyperfine interaction strength, suggesting that the hyperfine 
interaction may be intimately involved with the spin dynamics in this system. 

PACS numbers: 05.50.+q 75.30.Gw 75.40.Gb 75.50.-y 75.50.Lk 76.75. +i 



Ising models play a central role in our understanding of 
magnetic systems and their phase transitions. Their im- 
portance stems from their simplicity with respect to other 
models, and in the fact that they reproduce many ob- 
served physical phenomena (e.g. glassiness and quantum 
phase transitions). For T<2K and x<l, LiHo x Yi_ x F4 is 
thought to be a physical realization of the random trans- 
verse field Ising model with dipolar-magnetic interactions 
(plus a smaller nearest neighbor antiferromagnetic ex- 
change interaction) pHH. For 0.25<x<l the ground state 
of the system is a ferromagnet 0,0, Q; and for x<0.25 
enough randomness is introduced in the system such that 
long range ferromagnetic order is destroyed [1, 0, • ^ 
is natural to expect that in this last diluted regime, the 
long ranged dipolar interaction (which can be antifer- 
romagnetic for many bonds) together with the quenched 
chemical disorder produce a spin glass ground state. Sur- 
prisingly, the nature of the ground state for x<0.25 has 
been the topic of a heated debate at both the experimen- 
tal and theoretical levels. 

Experimentally it was found that for 0.1<x<0.25 the 
non-linear AC susceptibility (^3) peaks as a function of 
temperature, and that this peak gets rounded upon the 
application of an external magnetic field perpendicular 
to the Ising axis [7], |8(. These measurements were in- 
terpreted by the authors as a transition to a low tem- 
perature spin glass state. This interpretation was sup- 
ported by a numerical calculation 0, and the rounding 
of the peak in the presence of an external field was pro- 
posed to be a consequence of field induced random fields 
Other researchers though, pointed out that a crit- 
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ical analysis of \3 using only data above the peak (in 
equilibrium) indicates that there is no transition into a 
spin glass state at any finite temperature and transverse 
field IToL ITll . in agreement with previous numerical 
calculations jl2l . [l3L flJ ] . 

At a lower doping (x=0.045), one research group ob- 
served that the frequency response of the linear AC sus- 
ceptibility is narrower than that of a spin glass, and this 
was interpreted as a splitting of the system into clus- 
ters of spins which behave as single harmonic oscillators 
[l5l ] . In contrast, this narrowing was not observed in the 
measurements from another group which used a different 
sample with the same doping level. Instead, the tem- 
perature dependence of this xac data was shown to be 
compatible with that of a spin glass with a transition 
temperature lower than that achieved by the measure- 
ments [l6| . 

At this time there is no consensus on the nature of the 
ground state for x<0.25. One of the main reasons for this 
is the lack of experimental data with different probes, as 
available data is mostly on magnetic susceptibility and 
specific heat. With this in mind, we report in this let- 
ter muon spin relaxation (/zSR) [I7| measurements in a 
series of samples which span the whole diluted regime 
(x=0.25, 0.12, 0.08, 0.045 and 0.018). Our measurements 
do not present any of the signatures of canonical spin 
glass freezing; and an analysis of our data using a Kubo- 
Toyabe model [HI reveals that the Ho 3+ ions slow down 
with decreasing temperature down to «300mK. Below 
this point and down to wI5mK, the fluctuation rate of 
the Ho 3+ ions is observed to be temperature independent 
[igj | . Our analysis also shows that the low temperature 
fluctuation rate of the Ho 3+ magnetic moments increases 
with x up to 0.08 after which it levels off. 

Our samples are single crystals purchased from TY- 
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FIG. 1: Analysis of the low LF data using a power-exponential 
fitting function. The top panel shows the relaxation rate of 
the signal A, and the lower one the exponent j3. Color online. 



DEX J.S.Co. (St. Petersburg). Pieces from the main 
crystals were placed in the sample holder in such a way 
that the externally applied magnetic field was perpen- 
dicular to the Ising axis. The fiSK measurements were 
performed at the M15 and M20 beam lines of TRIUMF 
(Canada) in the Longitudinal Field (LF) configuration. 
In this configuration the initial muon spin direction is 
along the external magnetic field and therefore perpen- 
dicular to the Ising axis. For the measurements at M15, 
the samples were mounted on a silver sample holder us- 
ing "Apiezon N" grease for thermal contact. This holder 
was then screwed to the mixing chamber of a dilution re- 
frigerator. In this device the temperature of the samples 
was typically varied between 15mK and 3K, while the 
external field was scanned up to 0.2T. In the M20 beam 
line the temperature was controlled with a helium flow 
cryostat in the range from 1.8K to 100K, and the samples 
were mounted using a low background sample holder. 

Upon cooling from T«20K, the relaxation of the sig- 
nal increases due to slowing down of the magnetic Ho +3 
ions into the /iSR time window [H, H3] • The increase of 
the relaxation upon cooling is monotonic down to base 
temperature. Since a spin glass ground state is expected 
to be observed at these dilution levels, we analyzed our 
low LF data using a power-exponential fitting function: 
cxp(— (At)' 3 ). This function has been successfully used to 
study /iSR lincshapes of disordered spin systems, includ- 
ing spin glasses above T g [2l|. 

The fit values for A are shown in Figure [T] We will 
argue later that upon cooling, the increase in A around 
200mK (most noticeable in the x=0.25 data) is produced 
by a further slowing down of the magnetic moments of 
the system. This figure also shows the fit values for the 
power p. It can be seen that upon cooling, this param- 
eter has a minimum at around «10K, which is associ- 
ated with the slowing down of the very fast fluctuat- 
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TABLE I: Estimated freezing temperature (T m f), position 
of the peak of the specific heat (T pea k) [13], characteristic 
strength of the internal field (a or A), and low temperature 
fluctuation rate of Ho 3+ ions (fo), for each of the studied 
samples 



ing Ho moments into the /iSR time window [20| . Below 
this temperature, grows monotonically and it stabi- 
lizes at «0.85 for the x=0.12, 0.08 and 0.045 systems, 
and at «1.5 for x=0.25. If a glassy behavior was to be 
observed, this parameter should monotonically decrease 
upon cooling and reach a minimum of 1/3 just above the 
freezing temperature T g [2l| . An upper estimate for T g 
in LiHo x Yi_ x F4 can be obtained using the mean field 
expression: T m f « T g « xT c (T c is the critical temper- 
ature of the ferromagnetic x=l system which is 1.54K). 
These estimated temperatures are shown in Table HI to- 
gether with the temperature at which the magnetic spe- 
cific heat peaks [22| • Figure Q] shows that none of our 
samples presents a minimum of (3 at around T m f or T pea k- 
This fact together with the observation that no change in 
the shape of the signals is observed upon further cooling 
from «100mK, indicates that these systems do not have 
a canonical spin glass ground state. 

The low temperature data from the x=0.018 system 
could not be properly fit by a power-exponential since the 
/iSR signal develops a shoulder at low temperatures. This 
type of behavior indicates that fluctuations of the mag- 
netic moments in the x=0.018 system are slow. We de- 
cided to further analyze our data using Dynamical Kubo- 
Toyabe (DKT) polarization functions [1J|. The value of 
/? for the x=0.25 system at low temperature is approxi- 
mately 1.5 (Figure [T|) which indicates that, magnetically 
speaking, the system rather dense. Then, we used Gaus- 
sian DKT functions to fit this data. On the other hand, 
the low temperature values of (3 for the other systems are 
~0.85, indicating that these systems are magnetically di- 
luted. Then, we used Lorentzian DKT functions to fit 
the the data from the x=0.12, 0.08, 0.045 and 0.018 sys- 
tems. The fitting parameters of the Lorentzian (Gaus- 
sian) model are a (A) and v. a and A represent the 
characteristic size of the internal magnetic field at the 
muon site, while v is the fluctuation rate of the internal 
field (or the inverse correlation time of the local magnetic 
field, that is: (B(0)B(t)) cx exp(-vt)). We should men- 
tion that we attempted to analyze our data using other 
microscopic models (such as the the spin glass function 
in Reference (23l| ) b ut none them produced sensible (or 
physical) results |24| . 

The fit of the low LF data using the DKT functions 
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FIG. 2: Fitted values of v divided by its low temperature 
average value (vo) as a function of temperature and doping 
level. Color online. 

were performed by fixing a (A for the x=0.25 system) 
to the value found at base temperature, and then letting 
only the fluctuation rate v vary as a function of temper- 
ature. The fits with the DKT produced sensible results 
at the qualitative level [24[ . The values of a (A for the 
x=0.25 system) are shown in Table |U These values were 
found to roughly follow the -y/x trend expected for low 
values of x [24| . 

Figure [2] shows the value of v/vo as a function of tem- 
perature for all our samples, where vq is the average value 
of v in the low temperature range where it is constant. 
It can be seen that, as the temperature is lowered the 
fluctuation rate of the ions decreases; and below a tem- 
perature T* this fluctuation rate becomes temperature 
independent. T* does not seem to have a clear depen- 
dence on x, and has a typical value of 300mK. We note 
that this temperature is approximately the same size as 
the hyperfine interaction energy scale of the Ho 3+ ions 
(w200mK) [H. 

The fit values of v$ are shown in Table fl] It can be 
seen that Uq increases with Ho concentration until ap- 
proximately 0.08, above which it levels off. It is inter- 
esting to notice that the point of the x-T phase diagram 
where the characteristic strength of the dipolar interac- 
tion is smaller than that of the hyperfine interaction [26| 
is 0.13, which is close to 0.08, the point where we ob- 
served vq to flatten. 

The suitability of the DKT model to describe the low 
temperature behavior of LiHo x Yi_ x F4 was tested with 
the high LF data (see Figure [3]). Using the low LF 
fit parameters, we only increased the field in the DKT 
model and the resulting functions were observed to fol- 
low satisfactorily the experimental data up to «0.2T for 
the x=0.25, 0.12, 0.08 and 0.045 samples. The LF scan 
in Figure [3] clearly shows that the system is dynamic at 
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FIG. 3: Corrected asymmetry for LiHoo.12 Y0.88F4 at 12mK 
and different LF. Zero field - black squares, 0.01T - red circles, 
0.05T - blue triangles, 0.1T - green vertical squares, and 0.2T 
- violet stars. Lhe continuous lines are fits to Lorentzian DKT 
functions (a=12pis _1 , v=20[is~ 1 ). The partial decoupling for 
LF > 0.01T is due to a temperature-independent background 
present in the zero field data. The violet dashed line is the 
expected DKT line shape in 0.2T if the system was quasi- 
static (a=12/is -1 and ^=1.2/^s -1 ). Color online. 



low temperature. If the magnetic environment was static 
(frozen), it is expected that a LF of 0.2T would decou- 
ple the /xSR signal by about 85% (corresponding to an 
asymmetry of 0.14 in the figure). Instead, at this field 
the signal relaxes much below this point, evidencing the 
dynamic environment. Further more, the dashed line in 
the figure shows that the system is not even quasi-static. 
The fact that the /iSR signals for LF<0.2T are satisfac- 
torily described by the DKT model with a fixed value of 
Vq, indicate that the external field does not have a big ef- 
fect on this parameter (as expected since the Ising levels 
are not coupled for LF<2T Q). 

The experimental lineshape at intermediate LF 
(ss200G) in the x=0.018 system was not followed as 
closely as that for the denser systems. We believe that 
this is due to the internal field distribution of the system 
being stretched along the Ising axis [24]] > instead of the 
spherical one that the DKT function assumes. This de- 
formation has a big effect on the way the static signal is 
decoupled, and therefore its effects are most appreciable 
at this doping due to its low value of Vq (compared to the 
parameter a, see Table [XJ> . 

It should be remarked that our data shows the same 
qualitative behavior for all the samples over the explored 
temperature-field space. This is in contrast with the 
claim that the x=0.045 system has aphysically different 
ground state from the x=0.167 one Q. We should note 
though that specific heat and AC susceptibility measure- 
ments of our x=0.045 sample fltl I22I are different from 
those reported in References [5j, [l5|, |27| , opening the door 
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to the possibility that the the difference lies at the sample 
level. More measurements with other samples is needed 
to clarify this point. 

It is interesting to compare our results with those from 
other experimental techniques. Non-linear xac measure- 
ments in the x=0.198 and 0.167 systems show a peak at 
«140mK and «130mK respectively Q- These temper- 
atures, as well as those at which the magnetic specific 
heat peaks (Table HJ, are in the range where we observe 
the onset of temperature independent fluctuations, be- 
tween lOOmK and 600mK. Then, it is possible that the 
spin glass behavior of LiHo x Yi_ x F4 is associated with 
the slowing of the magnetic moments down to «300mK; 
but note that the temperature independent fluctuations 
at low temperature indicate that the freezing is not com- 
pleted (Figure [5]). This is in agreement with the absence 
of a spin glass transition obtained from a critical analysis 
of other non-linear %ac measurements 0, EH • 

It is natural to expect a spin glass ground state in 
LiHo x Yi_ x F 4 , since it posses the required characteristics: 
frustration, introduced by the dipolar interactions, and 
quenched disorder from the random dilution. As shown 
before, the ground state of the system is dynamic, and 
therefore is not a classical spin glass. In this sense, the 
observation of a dynamic ground state in LiHo x Yi_ x F4 
is as surprising as the observation of spin glassiness in 
the pyrochlore Tb2Mo2C>7 |28l 12911 fa similar compound, 
Y2M02O7, is also a spin glass |29l. l30l ]. but there is re- 
cent evidence for atomic position disorder in it [3lL I32I]). 
As the diluted dipolar Ising model does seem to have 
a spin glass ground state 0, it is tempting to believe 
that the hyperfine interaction can be responsible for the 
low temperature properties. The importance of this term 
for LiHo Y Yi _ x F4 has already been pointed out in Refer- 
ences [ILl33j: and we found that its energy scale coincides 
with the onset of the temperature-independent dynamic 
ground state. Also, as mentioned before, the hyperfine 



energy scale might also appear in the dependence of vq 
with dilution. Nevertheless, it is not clear how can the 
hyperfine term be responsible for the observed dynami- 
cal behavior, since this term has the effect of preventing 
fluctuations between the Ising levels instead of promoting 
them [l|. 

In summary, our /iSR measurements show that the 
ground state of LiHo x Yi_ x F4 for x<0.25 is not that of a 
canonical spin glass. We have observed that the low tem- 
perature state of the system is dynamical, and can be 
described satisfactorily by the stochastic Kubo-Toyabe 
model (which assumes a single correlation timescale). 
Using this model we have determined that the fluctu- 
ation rate of the Ho magnetic moments decreases as the 
temperature is lowered until «300mK. Below this tem- 
perature the fluctuation rate is temperature independent 
down to wl3mK. The /^SR signals from all our samples 
exhibit the same qualitative behavior as a function of 
temperature and doping which stands in contrast to the 
observation of an additional "anti-glass" phase inferred 
from some xac measurements. The hyperfine energy 
scale apparently manifests itself in the measured dynam- 
ical properties, suggesting that the hyperfine interaction 
is involved in the interesting dynamic behavior of this 
system. 
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